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ABSTRACT 

It has been recently reported that the electrical charge in a 
semiconductive carbon nanotube is not evenly distributed, but 
rather it is divided into charge “islands” [1,2]. This paper 
links the aforementioned phenomenon to tunneling and 
provides further insight into the higher rate of tunneling proc- 
esses, which makes tunneling devices attractive [3]. This pa- 
per also provides a basis for calculating the charge profile 
over the length of the tube so that nanoscale devices’ 
conductive properties may be fully exploited. 


1 INTRODUCTION 

Though it sounds futuristic, the idea of the warp tunnel is 
straightforward and sound, and it is essential for subsequent 
calculations in this paper. The concept of the warp tunnel is 
somewhat similar to Newton's Cradle. The warp tunnel is 
introduced in Figure 1. When the black mass in Figure 1 is 
inserted into the warp tunnel, its representative exits almost 
immediately from the opposite end of the tunnel. If the warp 
tunnel is full, the effective distance traveled by the black 
mass in “its journey across the tube” is on the order of its 
diameter. The warp tunnel's quantitative properties are de- 
rived here, based on Figure 2. m is the mass of the incoming 
or outgoing particle, and d its diameter. M is the mass of the 
average particle in the tube involved in the momentum trans- 
fer; D is the diameter of such a particle. The total mass of 
particles inside the tube is M (r/D). g' is the acceleration due 
to an arbitrary force F acting on m and the particles inside the 
tube. So, F = M (r/D) g\ assuming m « M (r/D). The equa- 
tions of uniform acceleration give d = g’ t 2 /2, t = (2 d/g’) 0 5 . 
Pretending that the outgoing particle is identical to the 
incoming particle, we can define m’s velocity across the tun- 
nel, v T , as 

v x = (r + d)/t = (r + d) (F/(2 d M (r/D))) 0 5 . 

Instead of F acting on a particle, we now assume it has a 
corresponding initial velocity Vj, then “it” will effectively 
traverse the tunnel with v T . Now we calculate an input 
velocity Vj equivalent to F based on energy considerations: 

mv^/2 = F d, so v 4 = (2 F d/m) 0 * 5 . 


Watp Tunnel 


As the particle is inserted into the warp tunnel, “it” 
simultaneously emerges from the other end, independent 
of the length of the tunnel. 


Axial Cross Section of Warp Tunnel 


The nearly simultaneous emergence of "the" particle as it is 
inserted is illusory. Actually, the outgoing particle emerges 
as a consequence of tunnel overflow and is not the selfsame 
particle inserted, but effectively represents it. 

Figure 1 : The warp tunnel. 



Figure 2: Auxiliary diagram for transmission velocity 
derivation. 


v-j/Vj = (r + d)/(2 d M (r/D)) 0 - 5 (m/(2 d)) 0 - 5 
= (r + d) (m D)°- 5 /(2 d)/(M r) 0 - 5 
a 0.5 (r/d) (m/M) 0 - 5 (D/r)° 5 , if r » d. 

An important conclusion of this analysis is that the incom- 
ing particle appears to traverse the tunnel with a velocity v T 
that is much greater than its input or initial velocity v { . The 
longer the tube is, the larger the mass, and the smaller the 
diameter or width of the inserted particle, the greater the 
speedup. However, due to conservation of energy, v 0 < v*. In 
this paper, we assume that a warp tunnel spans a potential 
barrier. This analysis then provides insight into the higher 
rate of tunneling processes. 

We derive inspiration from the optical tunneling scenario in 
Figure 3(a). When the upper right prism is removed, a detec- 
tor to the right of the lower left prism and in-line with the 
incident beam detects nothing. However, when the upper 
right prism is brought into place, that same detector now de- 
tects. Feynman discusses this phenomenon in ref. [4]. Though 
no translating photons are in the gap, oscillating electrons on 
the lower left prism create an oscillating field, which shakes 
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electrons on the upper right prism, thereby producing the 
light captured by the detector. 

We propose something similar for Figure 3(b). Electrons 
incident on the potential barrier do not actually enter the bar- 
rier, but rather produce photons free to enter the barrier in- 
stead. Through this action, electrons absorb energetic photons 
at the opposite end of the barrier. The energized electrons 
continue the motion once possessed by electrons on the other 
side of the barrier. In general, it is not the selfsame photon 
that enters the prism on the left that is detected beyond the 
prism on the right. Likewise, it is not the selfsame electron 
that impinges on the barrier in Figure 3(b) that emerges on 
the other side. 

We regard the neutral gaps between charge islands as being 
spanned by a warp tunnel. These warp tunnels are full of 
photons, so full that the injection of a single photon at one 
end of the tunnel causes the ejection of another at the other 
end near another charge island. The ejected photon is 
absorbed by an electron; maintaining, by this mechanism, a 
discontinuous charge distribution. 



2 SPACING CALCULATION 

Since charge is segregated into islands, to model this 
structure, an appropriate circuit seems to be that of a series 
RC circuit. A capacitor plate represents a charge island; the 
spacing between the plates of the capacitor represents the 
neutral gaps between adjacent charge islands. The capacitors 
are modeled as having equal capacitance while R represents 
the resistance of everything in the circuit except the nanotube, 
since the nanotube is essentially treated as a capacitive 
element in this circuit (Figures 4 and 5). Ordinarily, we 
would not expect a capacitor to conduct a DC current, but if 
the plate spacing is sufficiently close, a tunneling current will 
flow in the same way that electrons tunnel through thin 
barriers. When tunneling is viewed in this new light, a 
tunneling current can readily transmit through 40 nm barriers. 

Consider the warp tunnel in the top of Figure 1. Ignoring its 
coloring, let the particle entering the tunnel be a photon 
emitted by an electron, not displayed, to its left. Assume that 
the warp tunnel spans a barrier of width a. The particle 
emerging from the right side of the warp tunnel may be a 
photon that will be absorbed by an electron, not displayed, on 
the right side of the barrier. The outgoing photon thereby 
imparts that electron with the motion originally possessed by 


the other electron on the left side of the barrier. In these cal- 
culations, it is assumed that the photon has a width, denoted 
“It,” directly proportional to its energy or frequency. The 
proportionality constant a is derived below. 

Assuming the warp tunnel of length a in Figure 1 is full, 
the entering photon need only move into the tunnel an 
amount It when a representative photon will move out of the 
tunnel. The electron, which originally emitted the photon, has 
effectively bridged the distance across the barrier by “travel- 
ing” only a distance \£ « a. The effective length of the 
barrier is only It(f), the width or thickness of the atom of light 
of frequency f perpetuating the motion. 

Figures 4 and 5, which depict a DC RC circuit, motivate 
the derivation of tunneling current formulae. While ordinarily 
only a transient current flows in response to the application of 
a constant DC voltage across the series RC couple, given a 
certain dielectric between the capacitor plates, Giaever 
showed in ref. [5] that if do is small enough, a relatively con- 
stant tunneling current will flow, even for a DC power sup- 
ply. We can obtain this result from the classical analysis of a 
DC RC circuit provided the capacitor plate separation is d(t), 
a certain decreasing function of time, and where this d is a 
distinctly different variable from prior usages. From a 
mathematical perspective, the only d(t) that allows for a 

R’ - 37 MD 



Figure 4: The nanotube is modeled as a resistor and 
capacitor in parallel. However, the resistance R' of the 
couple is very high so that it may treated as an open 
circuit. R is the internal resistance of the circuit. It 
essentially includes resistance of the leads. R will be 
an essential parameter. 


Charge Islands 



semi-conducting 

dielectric. 


Regular cap. 
with dielectric; 



If d is above the nanoscale, the 
circuit current is transient, 
exponentially decaying. If d is on 
the nanoscale, a constant tunneling 
current exists. The classical ODE 
for an RC circuit may be used to 
calculate the tunneling current if 
the plate spacing d exponentially 
decreases with time. Equivalently, 
d may be fixed if we assume the 
density of transit lengths, L, goes 
as 1/L. 


Figure 5: Change islands modeled as charge on capacitor 
plates. 
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constant current in classical DC RC circuits is d(t) = d 0 e -s3 
where s3 is an undetermined constant. The capacitance C of 
the capacitor increases according to C = s/(dQ e _s3 r ) = s2 e s3 \ 
where s and s2 are constants easily determined by the defini- 
tion of capacitance. With that d(t), we obtain a constant tun- 
neling current through the capacitor, I c = (v/R) e _1/ ^ s2 s3) , 
which decreases exponentially with increasing initial barrier 
width d 0 . d(t) corresponds to a probability density of effective 

transit lengths L inside the plates (the originally fixed plates) 
that varies as 1/(L s3), a distribution that obviously favors 
shorter effective lengths. We assume, based on the model of a 
full warp tunnel, that the effective length of the barrier is 
about the width, It, of a photon, and that the more energetic or 
the higher the frequency of a photon, the thicker it is (It = 
a oj. For verification, the thermodynamic Planck Distribu- 
tion Law for the distribution of radiation of various frequen- 
cies inside a cavity may then give us the density of effective 
lengths inside the barrier. Using the approximation e x « 1 + x 
for small x, the Planck Distribution [6] simplifies to: 

df = l/(h to/(2 7 T t)), where h is Planck's constant and r is 
temperature. 

So the density goes as 1/to = l/[(lt/a)] x 1/L. This result is 
highly convincing, almost as much as the numbers that will 
be derived from it. Since the required distribution function for 
d(t) is actually df = 1/(L s3) [m _1 ], if L = It = ato, thermody- 
namics gives 

df = df = l/(h co/(2 7i i)> = l/(a co s3) 
yielding s3 = h/(2 % a x) [s' 1 ], where t is in Joules. 

To find a, now we derive an expression for the width of a 
photon as a function of its frequency. To derive an approxi- 
mation for the size of a photon, we assume that its density is 
similar to that of a neutron, which essentially has the maxi- 
mum known finite density in the observ able universe. 

We regard photons originating from both protons and elec- 
trons as having the density of a neutron. Since the proton is 
larger than the electron, the other charge quantum, it seems 
reasonable to approximate the upper size limit of one dimen- 
sion of the photon by the radius of a proton [7]. 

As shown in Figure 6, if in its rest frame a mass is spheri- 
cal, it will appear to observers in a rest frame viewing per- 
pendicular to the uniform motion of the mass as having the 
profile of a convex lens or vertical ellipse, due to its length 
contraction in the direction of its motion. As the relative ve- 
locity between the perpendicular observers' frame and the 
frame of the mass increases, the width of the mass will de- 
crease greatly while its height remains unaffected. In the 
limit, as the relative velocity approaches the speed of light, 
the mass will nearly have the profile of a rectangle. If observ- 
ers are in a rest frame looking in a direction parallel to the 
motion of the mass as it approaches (a dead center view), the 
mass will appear circular. So from the combined views of 
those observing perpendicular and parallel to the motion of 
the mass, the mass has a shape approximated by a hockey 


- 000(1 — 

Figure 6: Snapshot view of a mass in a rest frame perpen- 
dicular to its velocity. Masses towards the right have a 
successively higher relative velocity with the rest frame. 
Images of masses towards the right are laterally com- 
pressed due to the relativistic length contraction. Right 
most mass is depicted as being encapsulated by a rect- 
angle, which approximates its profile quite well. 

puck or a short, fat cylinder, which has both a rectangular and 
circular profile. This mass, the subject of this discussion, is a 
photon. The width of the rectangular profile or the short 
height of the hockey puck or photon is referred to here as the 
light thickness or It. The mass, mp, of the photon is given by 
mp c 2 = h f, where f is here used to represent the frequency as 
opposed to v. The volume of the photon, vp, is then the vol- 
ume of a cylinder, vp = It k r p 2 = mp/d’ = h f/(c 2 d'). 

lt(f) = (h f/(c 2 d"))/(7r r p 2 ). 

= 1.4IxlO" 39 f [m]. 

This result says that less energetic photons are thinner. We 
convert f into angular frequency: f = co/(2 7t). 

It(co) = 2.24x 1(U*° to [m]. 
a = 2.24x lO^ms. 

For the temperature variable, we use room temperature, 
298 K times Bk, Boltzmann's constant, 

t = 298 Bk = 4. 1 lx 10~ 21 J. 

s3 = h/(2 7t a t) = 1.1 5x 10 26 s _1 . 

From ref. [8] we use the low-end dielectric constant for 
carbon powder, 2.5. We assume a diameter of 1 nm for the 
nanotube to obtain its circumference, circ = 10 -9 n. We ob- 
tain 2.45 A from refs. [9] and [10] as the parallel carbon bond 
separation or the contribution of each atom along the CNT's 
cross section towards the circumference (Figures 7 and 8). So 
the number of atoms, na, along the CNT's cross-section is 
na = circ/2.45x 1 (H° = 12.82. 

The partial area, pa, an atom offers to its cross-section is 
obtained from refs. [9] and [10] as a double trapezoidal area 
with bases represented by the carbon bond length and the 
distance from opposite carbon atoms, and double height rep- 
resented by the parallel carbon bond separation. So, 

pa = ((2.83 + 1.42) x 1 (T 10 ) (2.45x 10“ 10 )/2 [m 2 ] 

= 5.21xl(T 20 m 2 . 

The cross-sectional area, a, of the circular wall of the CNT 
is a = na pa = 6.68xl0” 19 m 2 . The permittivity of free space 
£q is 8.85x 10 -12 F/m. 

s = kd 8 q a = 1.48x KT 29 F m. 
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Figure 7: Section of single non-helical shell. 


Equilibrium Structure 
Average Diameter of SWNT's 
Distance from opposite 
Carbon Atoms (Line 1) 
Analogous Carbon Atom 
Separation (Line 2) 
Parallel Carbon Bond 
Separation (Line 3)) 
Carbon Bond Length (Line 4) 


1.2- 1. 4 nm [1] 

S v 2.83 A [I] 

/ 2.456 A [1] 

2.45 A [1] 
1.42 A [1,23] 


Figure 8: Distances used in calculations. 


We may simplify the calculations by working with tunnel- 
ing resistances as opposed to tunneling currents. The tunnel- 
ing resistance of a capacitor is approximately given through 
this analysis as R T = R e Iy ^ R s2 s3 ->. 

Now we solve the tunneling equation to Find resistances 
corresponding to the measured separations. If the resistances 
are reasonable, our model is as well. Careful experiments 
could be performed to check these resistances, which we 
hope would be an outcome of this paper. We equate the tun- 
neling resistance of the entire nanotube with no islands at 
t = 0 when the power supply is switched on to that which 
would result from the known number n of charge islands 
which are manifested at t > 0: 

Exp[do/(R s s3)] = n Exp[d(/(n R s s3)], 


where in the first case, since the model says the charge is- 
lands should be evenly spaced, n = 650/40. In the second 
case, n = 750/36. 

Solving for R yields the internal circuit resistance for the 
given number of charge islands. The charge separations for 
various R have been calculated and are presented. 


3 RESULTS 


Results are given in the form: (R in Ohms, charge island 
spacing in nm). For the 650 nm tube, (0.0001270, 37.7), 
(0.0001316, 42.6), and (0.0001292, 40). For the 750 nm tube, 
(0.0001370, 34.4), (0.0001408, 37.82), and (0.0001389, 36). 
These resistances seem in range. They are very low resis- 
tances that appear to be what one would expect for the inter- 
nal resistance of the circuit under consideration. The charge 
island spacing increases with R. Note that these results differ 
from those reported in ref. [11]. There we accidentally used 
the double trapezoidal height in calculating pa, the incorrect 


value of 9.08xl0 16 kg/m 3 for the neutron density, an incor- 
rect value for the proton radius, and also a value of 5 instead 
of 2.5 for the CNT dielectric constant. 

4 CONCLUSIONS 

This model allows one to conclude that the effective length 
of the barrier may be given by the width, lt(f), of a photon 
within the barrier. To find the probability density function, 
f(L), of the length of the barrier, we use the Planck distribu- 
tion function, which gives the thermal average number of 
photons in a single mode of frequency. For the temperature 
variable, we use room temperature. The majority of the 
frequencies are smaller than that of visible light; hence, 
f(L) x 1/L, where L is a function of frequency. 

We return now to our classical RC circuit model. When the 
power is turned on in a classical DC RC circuit, we obtain a 
transient solution: a current that exponentially decays in time. 
However, if the capacitor plate spacing is not fixed, but is 
exponentially decreasing with increasing time, we find that 
the current through the capacitor is essentially constant and 
inversely exponentially depends on the original plate spacing. 
The probability distribution function, g(L), of the length L 
between the capacitor plates, or the plate spacing, must give 
g(L) x 1/L to explain a plate spacing that exponentially 
decays with time. This remarkable and highly convincing 
result shows that harmony may be obtained between standard 
tunneling current formulae and that predicted by DC RC cir- 
cuit analysis by assuming that the tunneled current follows a 
discrete trajectory [II] across the barrier. 
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